arXiv: 1509.02627vl [cs.CG] 9 Sep 2015 


Rectilinear convex hull with minimum area* 

Carlos Alegria-Galicia ^ Tzolkin Gardufio^ Areli Rosas-Navarrete^ 
Carlos Seara ^ Jorge Urrutia ^ 

September 10, 2015 


Abstract 

Let P be a planar set of n points in general position. We consider the problem 
of computing an orientation of the plane for which the Rectilinear Convex Hull of P 
has minimum area. Bae et al. (Computational Geometry: Theory and Applications, 
Vol. 42, 2009) solved the problem in quadratic time and linear space. We describe an 
algorithm that reduces this time complexity to 0(nlogn). 


1 Introduction 

The interest in the rectilinear convex hull of a set of geometric objects arises from the study 
of orthogonal convexity nanHi [20], a non-traditional notion of convexity that restricts 
convex sets to those (known as ortho-convex) whose intersection with any line parallel to 
a coordinate axis is either empty, or connected. For geometric object sets in the plane, 
the rectilinear convex hull has been extensively studied since its formalization in the early 
eighties, and has found applications in several research fields including illumination [T|, 
polyhedra reconstruction [9], geometric search [22|, and VLSI circuit layout design p3] . 
to name a few. Researchers have also studied relations between rectilinear convex hulls 
of colored point sets [3|, and explored generalizations of orthogonal convexity [nmaiiei 
along with related computational results HE]- 

We consider here the rectilinear convex hull as an orientation-dependent enclosing shape 
of a finite set of points in the plane. As every convex set is ortho-convex but not vice 
versa, the rectilinear convex hull of a point set is always contained in the standard convex 
hull of the same point set (and therefore, in several other traditional enclosing shapes), 
so it is relevant in applications where the enclosing shape is required to have minimum 
area. Applications with such requirements can be found in pattern recognition in digital 
images: rotation dependent and minimum area enclosing shapes are commonly used in 

*An extended abstract appeared in the XIV Spanish Meeting on Computational Geometry [2]. 

Rosgrado en Ciencia e Ingenieria de la Computacion, Universidad Nacional Autonoma de Mexico, 
{alegria_c, garduno_t, areli}-@uxmcc2. iimas .unam.mx. 

^Departamento de Matematica Aplicada II, Universitat Politecnica de Catalunya, Spain, 
carlos.seara@upc.edu. Research supported by projects Gen. Cat. DGR 2014SGR46, MINECO 
MTM2012-30951/FEDER 

Rnstituto de Matematicas, Universidad Nacional Autonoma de Mexico, urrutia@matem.unam.mx. Re¬ 
search supported in part by MTM2006-03909 (Spain) and SEP-CONACYT of Mexico, Proyecto 80268. 


1 



form shape analysis mm, and different enclosing shapes are continuously researched to 
be used as separators for feature classification [211 EH- 

Consider two rays £±,£2 leaving a point x G and an angle (j) such that, after counter¬ 
clockwise rotating £i by (p around x, we obtain £2. The polygonal chain £i U £2 divides the 
plane in two regions called wedges. We say that both wedges have apex x and sizes cj) and 
2tt — (j), respectively. Let Xg and Yg denote the x- and y- axis counter-clockwise rotated by 
an angle 6. A 9-quadrant is a translation of one of the (ortho-convex) wedges that result 
from subtracting Xg and Yg from the plane. We refer to a 0-quadrant as first, second, third, 
and fourth, according to the corresponding quadrant in the rotated coordinate system. 

Let P be a set of n points in the plane containing no three collinear elements and, for 
a fixed value of 9, let Rg{P) denote the (ortho-convex) region obtained by removing from 
the plane all first 0-quadrants that are P-frec, i.e., that contain no elements of P. The 
Rectilinear Convex Hull of P with orientation 0 is the set 
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nng{p) = f| Pe+,.|(P) 

1=0 


of points in the plane contained in all connected supersets of P that are ortho-convex with 
respect to Xg and Yg 13 EH]- See Figure 



(a) 




Figure 1: The rectilinear convex hull of P for three different orientations. The 
dotted polygon indicates the boundary of the standard convex hull of P. Note 
that TZ'Hg{P) is not always connected (disconnecting regions are bounded by 
the dashed rectangles in (b) and (c)). 


We solve the problem of computing all the values of 0 G [0, 27r) for which TZ'Hg{P) has 
minimum area in 0(nlogn) time and 0{n) space. Our general approach is to perform an 
angular sweep: we first discretize the set {0 : 0 G [0,27r)} into a sequence of increasing 
angles {0o,..., 0m-i}, m = 0{n), and then simulate the rotation of the plane by traversing 
the sequence in increasing order. While 0 runs from 0 to 27r, each 9i is an angle where there 
is a change in the relevant features of TZ'Hg{P) and, during every interval {9i, 0(j+i) mod m)-, 
those features remain unchanged. To sweep the sequence of angles we use an algorithm 
that:i) computes the angle 0 G [0i,02) that minimizes the area of Tl'Hg{P) in O(nlogn) 
time, and ii) obtains the angle of optimal area for the subsequent intervals in logarithmic 
time. We based our work on techniques from Bae et al. |7], Avis et al. |6], and Dfaz-Bahez 
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et al. m- For almost a linear factor, our result improves the 0{n?) time solution from Bae 
et al. [7] to an optimal 0(nlogn) time and 0(n) space. 

Outline of the paper. In Section we show how to compute and maintain the relevant 
features of TZ%g{P) during a complete rotation of the plane. In Section]^ we extend this 
result to compute the value(s) of 9 for which Tn-L 9 {P) has minimum area. Finally, in 
Section]^ we discuss our concluding remarks. 

2 The unoriented rectilinear convex hull 

In this section we introduce definitions that are central to our result. We also describe an 
algorithm to compute and maintain the relevant features of TZ'Hg{P) while 9 goes from 0 
to 27r in 0(nlogn) time and 0(n) space. 

2.1 Preliminaries 

As with every bounded ortho-convex region, the boundary of Rg{P) is a staircase-shaped 
ortho-convex curve. By construction, the curve in this case is an alternating polygonal 
chain with interior angle whose initial and hnal links are halflines parallel to Yg and Xg, 
respectively. In counterclockwise direction, every left-turn vertex in this polygonal chain is 
a point in P that we call a vertex of TZ'Hg{P). The staircase of P at angle 9 is the directed 
segment of the boundary of Rg{P), that starts at the rightmost and ends at the topmost 
vertex with respect to Xg and Yg (Figure [2(a)| ). Lets denote with Vg{P) the set of vertices 
over the boundary of Rg{P), ordered by appearance while traversing the staircase at 9 in 
its corresponding direction. For a fixed value of 9, we denote with 
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Vg{P) = U Ve+i-iiP) 

j=0 

the sequence of vertices of TZRg{P) in counterclockwise order. 



Figure 2: (a) Construction of a staircase, (b) The four staircases bound¬ 


ing TZHg{P). Staircases corresponding to opposite extremal 0-quadrants are 
marked with the same color. 
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At any particular value of 0, a P-free 0-quadrant is extremal if its apex is a clockwise- 
turn vertex in the staircase at 9. We say that extremal 9- and {9 -|- tt)- quadrants are 
opposite to each other (Figure 2(b)| ). Note that, if the intersection of two opposite 9- 
quadrants is not empty, then TZ'Hg{P) is disconnected. In such case, we say that both 
quadrants overlap, and refer to their intersection as their overlapping region. See Figure 



Figure 3: An overlapping region (filled with dark blue) generated by the inter¬ 
section of two extremal 9- and (0 -f tt)- quadrants (filled with light blue). 


Let Og{P) be the set of overlapping regions in Tn-Lg{P). For our purposes no specific 
order is needed in Og{P), so the overlapping regions can be stored in a constant-time query 
structure, such as a hash table. Remember that every staircase is monotone with respect 
to both Xg and Yg, and note that only one pair of staircases corresponding to opposite 
extremal 0-quadrants can intersect at the same time. Therefore, for a fixed value of 0, 
there is at most 0(n) overlapping regions in Og{P) that can be computed from Vg{P) in 
linear time. As for any fixed value of 0 the set Vg{P) can be compnted in 0(nlogn) time 
and 0(n) space [l5], we get the following theorem. 

Theorem 1. For a fixed value of 9, the (ordered) vertex set Vg{P) and the (not ordered) 
overlapping regions set Og{P) can he computed in 0(nlogn) time and 0{n) space. 

Note that, for a particular value of 0, each vertex of TZ'Hg{P) is the apex of a 0-quadrant 
free of elements of P. Conversely, a point in the plane lying in the interior of TZ'Hg{P) is 
the apex of (four) 0-quadrants containing at least one point of P. In a complete rotation of 
the plane (that is, by incrementing 0 from 0 to 27r), we thus recognize two types of events 
that modify Vg{P). An insertion (resp. deletion) event is a value of 0 where a point of 
P starts (resp. stops) being the apex of a P-free 0-quadrant. At every insertion (resp. 
deletion) event, a point of P joins (resp. leaves) V6i(P). 

Let p be an element of P. We say that a P-free wedge with apex at p is maximal, if it 
is bounded by rays each passing through another point from P. Note that every p that is 
the apex of a maximal wedge with size at least ^ is a vertex of TZT-Lg{P) during an angular 
interval in a complete rotation of the plane: the wedge can be transformed into an interval 
( 01 , 02 ) in such that 0i and 02 are respectively, the insertion and deletion events of p 
corresponding to (0i, 02) and, for any 0 G (0i, 02), P is the apex of a P-free 0-wedge; i.e., it 
is an element in Vg{P). See Figure]^ 

The vertex events table of P, is a list containing the set of angular intervals corre¬ 
sponding to all points in P that are apices of maximal wedges with size at least ^, ordered 
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Figure 4: Vertex events of the point p G P. The dotted lines denote the 
bounding rays of 0-quadrants: (a) p is the apex of a maximal wedge with size 
at least (b) p is the apex of a 0-quadrant containing at least one point of 
P. (c) Insertion event. |(d)| p belongs to Ve{P)- (e) Deletion event. |(f)| p is no 
longer a vertex of TZHeiP)- 


by appearance as 9 goes from 0 to 27r. As every p can be the apex of at most three of 
such wedges, the number of vertex events (and the angular intervals they delimit) is upper 
bounded by 0(n). 

Theorem 2 (Diaz-Bahez et al. [IH])- The vertex events table of P can be computed in 
0(nlogn) time and 0{n) space. Using this table, Ve{P) can he maintained over all 9 G 
[0,27r) in O(relogn) time and 0{n) space. 

In general terms, in the algorithm of Theorem the vertex intervals table is used as 
follows. For an initial value of 9, four vertical lines are located in the table at angles 9, 
0 -b f, 9 + TT, and 9 -|- as shown in Figure The intervals intersected by each line 
correspond to points in P that are apices of first, second, third, and fourth 0-quadrants 
free of elements of P and thus, are stored in the ordered sets Vq{P), Vq^k(P), Vg+Tt{P) and 
respectively. The rotation of the plane is then simulated by sweeping the table 
towards 27r using the four lines at the same time. A line over an interval endpoint marks 
a value of 0 when a vertex event occurs, so Vg{P) is updated accordingly. We refer the 
reader to Dfaz-Bahez et al. |T0] for more details on the sweeping process. 

During a complete rotation of the plane, we identify two types of events that modify 
Og{P). An overlap (resp. release) event is the value of 0 where two opposite 0-quadrants 
start (resp. stop) generating an overlapping region. At an overlap (resp. release) event. 
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Figure 5: Using the vertex events table to compute the vertices of 'R'He{P)- 


an overlapping region is added to (resp. deleted from) O 0 {P) (Figures and [^. As we 
mentioned before, this updates can be done in constant time as we do not require the 
overlapping regions in Og{P) to follow any particular order. 







Figure 6: The life of an overlapping region: (a) A pair of extremal 9- and (0+7r)- 
quadrants are about to overlap. m Overlap event. |(c)| The overlapping region 
belongs to Og{P). Ea Release event. Note this event matches the deletion 
event of the highlighted vertex. I (e) I The overlapping region no longer exists. 


To maintain TZ'}i 0 {P) through a complete rotation of the plane, we will use the sets 
V 0 {P) and Og{P). The set Ve(F’) will be maintained using the sweeping process from 
Theorem To maintain O 0 {P), during the rest of this section we will show how to 
compute an ordered list of overlapping events, that resembles the vertex events table. The 
idea is to sweep both tables at the same time, while updating Ve(P) and O 0 {P) at each 
event. As can be seen in Figure overlapping events do not necessarily match up with 
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vertex events and thus, both tables have to be computed independently. 





(a) 


(b) 


(c) 


Figure 7: A pair of overlapping events not matching up with any vertex event. 


2.2 Characterizing overlapping regions 

As we increment 6 from 0 to 27r, the apex of every extremal 0-quadrant traces at least 
one oriented circular arc. This arc is a segment of a semicircle whose diameter is the line 
segment joining the vertices of TZHeiP) supporting the particular quadrant. Let a be one 
of the arcs traced by an extremal 0-quadrant and ip the angular interval such that, for all 
0 G if, the apex of 0-quadrant belongs to a. We say that a is an extremal arc of P, and 
refer to ip as its tracing interval. See Figure 








Figure 8: Extremal arcs construction for two point sets. The “x” marks denote 
common endpoints between arcs. Tracing quadrants and their corresponding 
arcs are marked in blue. (a)|(c) An extremal arc sequence. (d)|[(f) A quadrant 
supported by the same pair of vertices can trace several extremal arcs. 


Consider two extremal arcs ai and 02 and their corresponding tracing intervals ipi and 
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(^ 2 - Assume the arcs are traced by a pair of opposite extremal 0-quadrants supported by 
the points a,b, c, d £ P as shown in FigureThe full overlap interval of ai and 02 , is the 
angular interval [9i , ^ 2 ] where Xq^ and are respectively, parallel to the line segments 
ac and bd. The facing interval of ai and 02 is the (possibly empty) angular interval given 
by ifi n {ip 2 + tt). Note that ai and 02 may admit overlapping quadrants only if their full 
overlap interval is proper {9i < 02)) and their facing interval is not empty. If this is the 
case, the extremal 0-quadrants overlap during at most one angle interval [7], that can be 
obtained by intersecting the full overlap interval with the facing interval of both arcs. The 
following lemma is straightforward. 

Lemma 1 (Overlap in terms of arcs). Two extremal arcs admit overlapping quadrants 
if, and only if, their full overlap interval is proper, and the intersection between their full 
overlap interval and their facing interval is not empty. 





Figure 9: Full overlap intervals: (a)|(b) A proper full overlap interval (0i < 02)- 
(c)|(d) A non proper full overlap interval (0i > 62 ). 


For a particular value of 0, a vertex of TZTi 0 {P) supports at most two extremal 0- 
quadrants, so each insertion (resp. deletion) event generates (resp. terminates) a constant 
number of arcs. The number of vertex events and extremal arcs thus share the same linear 




upper bound. As the number of overlapping intervals is linear in the number of points 
in P [7], by testing Lemma on every pair of arcs we can obtain in O(n^) time the set 
overlapping intervals, sorted by appearance while 9 goes from 0 to 27r. We devote our 
following discussions to lower down this brute-force complexity, by reducing the number 
of candidate pairs of arcs to 0(n). 

Let C7i(P) be the (standard) convex hull of P, {pi, ... ,ph} the set of points in P over 
the boundary of C'H{P) labeled in counterclockwise order, and {ei,..., Ch} the edge set of 
C'}i{P), where Cj = pipi+I for 1 < i < /i — 1 and = PhPi- We denote with Ai{P) the 
sequence of extremal arcs ordered by their corresponding direction, that were traced by 
extremal 0 -quadrants corresponding to every segment of a staircase connecting pi and pi+i 
at some 9 during a complete rotation of the plane. We further define the ordered sequence 
of extremal arcs A{P) = Ai{P) U A 2 {P) U • • • U Ah{P) and say that A{P) and Ai{P) are 
respectively, the arc chain of P, and the subchain of P corresponding to Ci (see Figure [Tol). 



Pi+i 


Figure 10: The arch-chain of P, and the subchain Ai{P) marked in blue. 


Note that Ai{P) is a continuous curve whose endpoints are the endpoints of e*. Every 
point a G Aj(P) is the apex of a P-free wedge Wa, bounded by rays passing through some 
point in P. This wedge have size at least f if a belongs to P m and exactly | otherwise, 
as it is an extremal 0 -quadrant. 

Lemma 2 (Avis et al. [6j). Let I be the line orthogonal to ei passing through a. The rays 
of Wa intersect e* in opposite sides of i (^Figure [17]). 

Corollary 1. No point b A ct in- Ai{P) belongs to Wa- 

Proof. Suppose there is a point 6 7 ^ a of Ai{P) contained in Wa- If b belongs to P, then 
Wa is not P-free, and we are done. If this is not the case, then b is the apex of a P-free 
0-quadrant qb, with bounding rays passing through some point of P. By Lemma at 
least one of the rays of Qb is completely contained in Wa- Therefore, at least one of the 
points of P supporting the rays of qb is also contained in Wa and then, Wa is not P-free (a 
contradiction). □ 

Lemma 3 (Monotony). Any line orthogonal to Ci intersects Ai[P) at most once. 

Proof. Suppose there is a line orthogonal to e* that intersects Ai{P) more than once. Let 
a be the intersection point farthest from Cj, and Wa its corresponding P-free wedge. By 
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Figure 11: The subchain corresponding to e^. 


Lemma[^ any intersection point a' ^ a\s in the interior of Wa- In contradiction, Corollary[^ 
states that a' does not belong to tCa- □ 


Consider the orthogonal projections of the points in Ai{P) over Cj. We denote with -<i 
the following relation between any a,b £ Ai{P): a 6 if the projection of a is closer to 
Pi than the projections of b. We use a Pi b to denote that either a Pi b, or both points 
have the same projections. By the monotony lemma, Pi is a total order in the elements 
of Ai{P) and thus, a is equal to b if their projections coincide. Let (pj,... ,pIh) be the 
sequence of points in P belonging to Ai{P), where pi = p\ Pi P2 Pi • • • Pi Phi — Pi+i- The 
link between pt and is the set of points 


= {a £ Ai{P)\ p] Pi a Pi p*+i} 


that form a segment of Ai{P) whose endpoints are p*- and Pj+i- See Figure 
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Figure 12: The link Cj- C Ai{P) marked in blue. Pointed line segments denote 
the vertical projections of the points in P over e^. 


Note that all the segments of a staircase that traced Ai{P) during a complete rotation 
of the plane, are contained in the P-free region bounded by the closed curve Ai{P) U 
Ci- Therefore, each pair of extremal 0-quadrants generating an overlapping region, trace 
arcs that belong to links that intersect (please refer again to Figure 10). From the O(n^) 
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total pairs of extremal arcs, we thus discriminate those belonging to intersecting links as 
candidates to admit overlapping 0-quadrants. 


2.3 A linear number of intersections between links 

In order to compute all pairs of intersecting links in P, we will first obtain a bound on 
the number of their intersection points. After several auxiliary results, in this section we 
present a proof for the following theorem. 

Theorem 3. There are 0{n) intersection points between the links eontained in A{P). 


Upper bounding intersecting links. Our general approach to prove Theorem is to 
enclose links and their corresponding subchains in disks. We then apply combinatorial 
properties of the resulting disk arrangement from Avis et al. [6|, to upper bound the 
intersecting pairs of links to a linear number. We will start by proving the cornerstone of 
our approach. 


Lemma 4. Let a, b, c be three pairwise different points in iij such that a <ib <i c. Then, 
I < Zabc < vr. 


Proof. Note that b does not belong to P and thus, it is the apex of a P-free 0-quadrant 
bounded by rays passing through some point in P. By Corollary neither a nor c are in 
the interior of this quadrant and therefore, Za6c > |. Let p\ and be the points of P 
supporting the rays of the quadrant at b. Without loss of generality, assume that pi -<i p)^ 
as shown in Figure 13 By Lemma it is easy to see that 


Pi Zi p] Pi a <ib <i c Pi pj^i Pi 


The points a and c are vertices of P-free wedges with size at least ^, so they are contained 
in the disk D with the line segment as diameter. As b lies over the boundary of D, 
we conclude that Pabc < vr. □ 


b 



Figure 13: Construction for Lemma 


0 
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The edge-disk corresponding to Cj is the disk Di having e* as diameter. The link-disk 
corresponding to iij is the disk Dij with the line segment as diameter. We denote 

with Tj and rij the corresponding radius of Di and Dij and, for the sake of simplicity, 
we abuse terminology by saying that rjj is the radius of Note that by Lemma 
Ai{P) C A and lij C A,j- 

Lemma 5. A link-disk Dij is intersected by at most 28 edge-disks with radius greater than 
or equal to rij. 

Proof. A set P formed by m disks is a k-set, 0 < /c < m, if no point in the plane belongs 
to more than k disks in P. From Avis et al. [ 6 ] (specifically. Section 2.1) we know that, in 
a k-set, the disk with the smallest radius is intersected by at most Ik disks. Moreover, it 
is also shown that all the edge-disks induced by P form a 3-set. 

Consider now the set P{P) = {A : fg > r-jj}. As P{P) is a 3-set, the set P'{P) = 
'D(P)U{A,j} is a 4-set where Aj is the smaller disk. The link-disk Aj' is thus intersected 
by at most 28 edge-disks m.P'{P). □ 


Lemma 6. Let Dij be a link-disk intersected by some edge-disk Dg. At most five links in 
Dg with radius greater than or equal to rij intersect Dij. 


Proof. Let R be the strip-shaped region bounded by the lines that orthogonally project 


Dij over Cg 


From Lemma the links in Dg that intersect Dij must have at least one 


endpoint contained in R (see Figure 14). 



Figure 14: Links intersecting Di j must have at least one endpoint in R. 

Lets denote with df the minimum Euclidean distance between pf and Cg, and say that 
pI is a peak (resp. a valley) in As{P), if df is greater than (resp. less than) both df_i and 


12 




Lets suppose there are three points contained in R, such that pf is a 

peak in As{P). Without loss of generality, assume that df_i > df_^_l and, for simplicity, 
that all the links in Dg with at least one endpoint contained in R, have radius greater than 
or equal to rjj. As rg^t-i > the point pf does not belong to the disk marked as D in 
Figure 15(a)). Specifically, D is contained in R, has radius rjj, and its diameter is a line 
segment parallel to e* that passes through pl_i. 



Figure 15: (a) There are no peaks in R. (b) 



Note that p| cannot be a point in As{P) as in contradiction with 

Lemma 1^ Therefore, the region R contains no peaks and thus, at most one valley. As at 
most two link endpoints from As{P) can be antipodal in Dij, it is easy to see that at most 
five links in Dg intersect Dij (Figure 15(b)). □ 


Lemma 7. There are 0{n) pairs of intersecting links in A{P). 


Proof. Let C be the list of all the links induced by P ordered by increasing radius. From 
Lemmas[^and[^ the first link in L is intersected by at most a constant number (specifically, 
28 X 5 = 140) of the remaining links in C. By removing this link from T, we get that the 
next link in the list is also intersected by at most constant number of links. As there is a 
linear number of extremal arcs and each arc belongs to a single link, there is also a linear 
number of elements in C. Therefore, by recursively removing the link with smallest radius 
from T, the total number of intersecting pairs adds up to 0(n). □ 


The proof of Theorem 

Proof. From Lemma we know that two links within the same subchain intersect only 
at one of their endpoints. We will prove next that links belonging to different subchains 
intersect at most twice. The theorem follows from this result and Lemma 0 
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Suppose there exist at least three intersection points a, b, c between two links iij and 
£s,t and, without loss of generality, assume that a <ib <i c. Observe that a <sb <s c\s the 
only possible order of the projections of these points with respect to e^. For any different 
order, some of these points cannot be contained in (.s,t because they form an angle either 
smaller than | (Figure 16(a)) or greater than vr (Figure 16(b)), contradicting LemmaHl 



Lets assume that a 6 c and denote with e; and Cm, the edges of C'H{P) intersected 

by the rays from b passing through a and c, respectively. Note that have to be such that 


m < s < i, or i < s < 1. See Figure 17(a) 



Suppose that m < s < i. Let i be the line perpendicular to passing through a, and 
Wa the wedge corresponding to a. Remember from previous discussions that Wa is P-free, 
have size at least is bounded by rays intersecting e^, and does not contain any other 
point from £ij (see Figure [T7(b)[ ). Note that c and Ps+i are in opposite sides of £ and are 
not contained in Wa- Therefore, | and Zacps+i < f • As by our hypothesis 
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a c Ps+i, we get from Lemmaj^that c is not contained in As{P). If i < s < Z, by a 
similar argument we can deduce that a cannot belong to As{P). □ 

2.4 Computing the overlapping intervals table 

Let us color the extremal arcs of P as follows: red if it belongs to a subchain corresponding 
to an edge in the upper chain of C4Z(P), and blue otherwise (see Figure [T^. Note that, 
regarding of the value of 6, any pair of extremal 0-quadrants intersecting an edge in the 
upper chain (or lower chain) of C'H{P), are not opposite to each other. The facing interval 
of any pair of monochromatic extremal arcs is thus empty and therefore, they do not admit 
overlapping 0-quadrants. 



Ui+ll 


Figure 18: The colored arc chain of P, and the labels assigned to a pair of sets 
of consecutive extremal arcs. 

Let (ai, 02 ,..., ao(n)) be the sequence of extremal arcs forming A{P), labeled by ap¬ 
pearance while traversing the arc-chain following the extremal arcs corresponding direction 
(see Figure We denote with iu,v the subsequence (a^, Om+i . •., ay-i, a„) of consecutive 
arcs in A{P) forming a link. Note that by construction (see again Figure]^, the extremal 
arcs in iu,v define a sequence au,v = {oiu, ■ ■ ■ of increasing angles where [ai,ai+i], 

u < i < V, is the tracing interval of a*. 

Consider iu,v and £s,t to be a pair of intersecting links. In 0{nu,v + ns,t) time, where 
nu,v = V — u + 1 and = t — s + 1, we can merge au,v with the sequence 

— ^s,t -h vr = {ug -|- vr,..., at+i + tt) 

into a sequence 7 = ...) of increasing angles, where the jth element ai. comes 

from au^v A u < ij < v, or from a'g^ if s < ij < t. Assume the smallest element in 7 is an 
angle coming from Let a^- and be to elements in 7 such that ai^ comes from au,v, 
ai^. comes from a'g^, ij is greater than ik, and there is no G 7 such that ij > ik' > ik 
and s < ik' < t. In other words, ai^, is the element in 7 immediately preceding Ui- that 
comes from a'g^. Note that 

n [0:4, 7^0; 

that is, the extremal arcs ai- and aj^. have a non empty facing interval. By performing 
a linear run in 7 , we can therefore discriminate from all the pairs of extremal arcs in 
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the intersecting links, those with non empty facing interval. In a second linear run, we 
can use Lemma to obtain the 0{nu,v + overlapping intervals of the arcs forming 
the intersecting links. Following a similar reasoning we can deal with the case where the 
smallest element in 7 comes from au,v We get the following lemma. 

Lemma 8. Let iij and be two intersecting links formed by Uij and extremal arcs, 
respectively. The extremal arcs in iij and ig^t generate 0{nij + ng^t) overlapping intervals, 
that can he computed in 0{nij + Ug^t) time and space. 

We now outline the algorithm to compute the overlapping intervals table. 

1. Compute the vertex intervals table. 

Use the algorithm from Theorem 

2. Compute the arc chain of P. 

Each extremal arc should be specified in terms of its tracing interval, and the points 
supporting the corresponding extremal 0-quadrant. The elements in A{P) should be 
ordered according to the arcs respective direction and grouped by links. These re¬ 
quirements can be fulfilled using the vertex intervals table, by extending the traversal 
algorithm from Diaz-Bahez et al. do] we mentioned in Theorem § 

(a) At each insertion event, at most two extremal arcs are generated, and at most 
one extremal arc is interrupted. Pointers should be set up from the interrupted 
to the newly created extremal arcs. If an extreme of a new arc is a point in P, 
a new link should be initialized with the respective extremal arc. 

(b) At each deletion event, at most one extremal arc is generated, and at most 
two extremal arcs are interrupted. One of the interrupted extremal arcs will be 
always ending at a point in P, so a link is completed. As before, pointers from 
the interrupted to the newly created extremal arcs should be set up. 

3. Color extremal arcs. 

Color each arc with either blue or red as discussed previously in this section. This can 
be done by traversing A{P) while assigning colors to each arc, according to whether 
the current edge e* belongs to either the upper, or the lower chain of CTL{P). 

4. Identify bichromatic intersecting links. 

Note that the largest possible extremal arc is a semicircle and thus, it should be 
partitioned in at most three segments to get a set of curves monotone with respect 
to an arbitrary direction. The arcs in A{P) can thus be transformed into a set 
A' (P) of curves monotone with respect to the same direction, having the same upper 
bound on the number of elements as A{P). We can therefore apply the Bentley and 
Ottmann plane sweep algorithm [ 8 | on A'{P) to compute the intersection points be¬ 
tween extremal arcs. From these intersection points, we discriminate those belonging 
to bichromatic pairs of arcs. Pointers to the links containing the involved extremal 
arcs should be setup, so we can obtain the set of all bichromatic intersecting links. 
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5. Compute the overlapping intervals table. 

Apply the algorithm from Lemmaj^on each pair of intersecting links. The overlapping 
intervals should be stored in a list, ordered by appearance while rotating the plane 
from 0 to 27r. Along with each overlapping interval, information to construct the 
corresponding overlapping region should also be stored. 

By Theorem]^ item requires 0(nlogn) time and 0{n) space. As there is a linear 
number of vertex events and a constant number of additional operations are performed 
at each event, item also consumes 0(n log n) time and 0{n) space. Item takes 0{n) 
time/space, as there is a linear number of extremal arcs in A{P). To compute the set 
A'{P) we require a linear run on A{P). By Theorem]^ the Bentley and Ottmann plane 
sweep algorithm processes A'{P) in 0(n logn) time and 0{n) space, so item|^also requires 
0(n log n) time and 0{n) space. Finally, from Lemmas and and the fact that there 
is a linear number of extremal arcs and each arc only belongs to one link, item requires 
0(n logn) time and 0{n) space. 

Theorem 4. A table containing the (linear number) of overlapping intervals induced by all 
points in P sorted by appearance while 6 goes from 0 to 27r, can be computed in 0(n log n) 
time and 0{n) space. 

2.5 Maintaining TZ'Hg{P) 

The overlapping and the vertex intervals tables can be merged in 0{n) time into a single 
table containing a linear number of intervals ordered by appearance while 0 goes from 0 
to 27r. The extremes of these intervals effectively discretize the set {9 : 9 € [0,27r)} into 
a sequence {0i,..., 0o(n)} of increasing angles. While 9 varies in [0, 27r), TZ'}iQ{P) remains 
unchanged during each interval (0j,0i+i), and each angle 9i is an insertion, a deletion, an 
overlap, or a release event. 

At any 9 G [0i,02), we compute the sets Ve{P) and Oo{P). We then simulate the 
rotation of the plane by sweeping the events table towards 27r, using four vertical lines 
at angles 9, 9 + ^, 9 + tt, and 9 + jTt, in a similar way as we described in Section ^ 
When a sweeping line reaches a vertex event, the set Ve{P) should be updated accordingly. 
Overlapping regions are detected by sweeping lines at 9 and 0 + f reaching an overlap 
event. The line at 9 detects overlapping regions from staircases at angles 9 and 9 + tt, 
and the line at 0 + | the ones belonging to staircases at 0 + | and 9 + |7r. At each 
overlapping event, the set Og{P) should be updated accordingly. If Og{P) is stored in a 
hash table (an ordering in the overlapping regions is not mandatory), these updates can be 
done in constant time. Considering previous discussions on this section and in particular, 
the complexities of Theorems 00 and and the fact that the events table has a linear 
number of elements, the maintenance procedure we have just described requires 0(n log n) 
time and 0(n) space. Note that this bound is tight as at any fixed 9, C'H{P) can be 
obtained from Vg{P) in linear time. 

Theorem 5. The sets Vg{P) and Og{P) can be computed and maintained while 9 goes 
from 0 to 27r, in 0 (nlogn) time and 0{n) space. 
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3 Minimum area 


In this section, we adapt characterizations from Bae et ah [7] to the idiosyncrasies of the 
algorithm of Theorem We then show how to extend this algorithm to compute the value 
of 9 for which the area of TZl-Le{P) is minimized in 0(nlogn) time and 0{n) space. 

Consider again that {0i,..., is the sequence of (vertex and overlapping) events, 

ordered by appearance while 9 goes from 0 to 27r. We express the area of TZ'Ho{P) for any 
9 G [9i,9i+i) as 

area(7^He(P)) = area('P(0)) — area(Aj(0)) + area(Dfc(0)), (1) 

j k 

where V{9) denotes the polygon obtained by connecting with line segments consecutive 
vertices over the staircases forming TZT-Lq{P), Aj{9) is the jth triangular region bounded 
by the union of the line segment and the segment of an extremal 0-quadrant both connecting 
a pair of consecutive vertices on a staircase, and □fc(0) is the kth overlapping region defined 
by the intersection of two opposite extremal 0-quadrants. See Figure [791 



Figure 19: Computing the area of Tl'Hg{P). The polygon P{9) is bounded by 
the dotted line. A triangle Aj(0) and an overlapping region □^(0) are filled 
with blue. 


Our general approach is to maintain Equation 0 during a complete rotation of the 
plane, by traversing the event sequence while adding and removing terms accordingly 
at each event. At the same time, we obtain the angle of minimum area of TZ7iQ{P) by 
computing the optimal angle of each interval [0j,0j+i), and replacing the global optimal 
angle only if the minimum area is improved. As we show next, at any particular value of 
0 we can evaluate Equation Q in linear time, and a constant number of terms need to be 
updated at each event, regarding its type. 


The polygon. At any fixed value of 0, the area of V{9) can be computed from Vo{P) in 
linear time. During an interval between events the polygon remains the same, so area(7^(0)) 
only changes at vertex events. These events can be handled in constant time: the area of 
a triangle needs to be subtracted (resp. added) to the previous value of area('P(0)) at an 


insertion event (resp. deletion event). See Figure 20 
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(a) 


(b) 


Figure 20: (a) The point p is about to become a vertex, (b) After the inser¬ 


tion event of p, the area of the triangle in blue needs to be subtracted from 
area('P(6i)). 


The triangular regions. A triangular region is defined by a pair of vertices in a staircase 
forming TZ7iQ{P), that are consecutive in the staircase’s corresponding direction. If we 
consider a staircase at angle 9 and two consecutive vertices Vj,Vj+i £ Ve(P), the area of 
Aj(0) is bounded by a line through Vj and Vj+i, a line through vj parallel to X^, and a 
line through Uj+i parallel to Yq. From Bae et al. [7], we know that the area of the jth 
triangular region can be expressed as 


area(Aj(0)) = cj cos 2(0 — 9i) -|- dj sin 2(0 — 0j), 


where Cj and dj are constant values depending on Vj,Vj+i, and 0*. The term area(Aj(0)) 
is only impacted by vertex events. At each event, it has to be modified a constant number 
of times. At an insertion event, the term for one triangle is removed and, as any vertex 
supports at most two extremal 0-quadrants, the terms of at most two triangles need to be 
added. The converse occurs for deletion events. See Figure 




Figure 21: Updating the term area(Aj(0)). At an insertion event, (a) 
most one triangle should be removed, and m at most two triangles should 
added. 


at 

be 


The overlapping regions. An overlapping region is a rectangle isothetic with respect 
to Xg and Yg. If we consider a pair of opposite extremal 0-quadrants, from Bae et al. [7] 
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we know the area of the kth overlapping region can be expressed as 

area(nfc(0)) = Ck + dk cos 2(0 - 0j) + sin 2(0 - 0*), 

where Ck,dk, and are constants depending on 0j and the vertices supporting the over¬ 
lapping quadrants. The term area(nfc(0)) is impacted by all types of events. Overlap 
events require the term of a single overlapping region to be added or deleted. Again, as 
any vertex supports at most two extremal 0-quadrants, at a vertex event the terms of a 
constant number of overlapping regions are added, or deleted. Please see again Figure 
for reference. 


Characterization. Before describing our algorithm, we consider some basic questions 
about the behavior of aTea,(TZ'He{P)). First of all, from Lemma 4 of Bae et al. there 
seems not possible to restrict the number of candidate angles of minimum area. Second, 
within an interval between events, aie.a(JZl-L 0 {P)) has a linear number of terms (the number 
of triangles and overlapping regions share the same upper bound as the elements in Vb{P)) 
and thus, we can reduce it to 

c + Dcos2{e-ei) + Es\ii2{e-ei) ( 2 ) 

in 0{n) time, where C, D and E are the constants resulting from adding up the constants 
of the terms in Equation Q. As Equation Q has a constant number of inflection points 
in [0, 27r), an additional constant number of operations allow us to obtain the local angle 
of minimum area within the specific angular interval. 


The search algorithm. We now outline the algorithm to compute the angle of minimum 
area. 


1. Compute the sequence of events using the algorithm from Theorem]^ 

2. At the first interval of the sequence, compute V 0 {P),Og{P), area{TZ'Hg{P)) as ex¬ 
pressed in Equation Q, and the angle 0m of minimum area. 

3. Traverse the sequence of events. At each event: 


(a) Update the sets V 0 {P) and O 0 {P). 

(b) Handle each event updating Equation Q by adding or subtracting terms as 
explained before in this section. 

(c) Compute the local angle of maximum area using the updated version of Equa¬ 
tion (§. Replace 0m only if the current minimum area is improved. 


Erom Theorem]^ item[^takes 0(n log n) time and 0{n) space. Itemj^requires 0(n log n) 
time and 0{n) space (Theorem [^, plus an additional linear run to obtain Equation Q 
and the initial value of 0m- As there is a linear number of events in total, items 3a to 
are repeated 0(n) times. Item consumes logarithmic time (see the sweeping process of 
Theorem]^ and, as we previously described, items 3b and [3^ take constant time. Einally, 
note that a trivial modification in this algorithm allow us to report the angles of minimum 
area, if there are more than one in a complete rotation of the plane. Erom this analysis we 
obtain our main result. 
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Theorem 6. Computing the value(s) of 9 € [0, 27r) for which TZ'Hg{P) has minimum area, 
requires 0 (nlogn) time and 0{n) space. 

4 Concluding remarks 

We used a preprocessing step created by extending a known 0(n log n) time algorithm, 
to maintain TZ'Hg(P) as 6 changes from 0 to 27r. Using this result, we described how the 
orientation of the plane in which TZT-Lg[P) has minimum area, can be obtained in 0(nlogn) 
time and 0{n) space. If there exist more than one optimal orientation, the algorithm finds 
them all with no additional complexity in time nor space. 0 ur result improves the quadratic 
time solution presented by Bae et al. [7] to an optimal 0(nlogn) time and 0(n) space. 

Notice that our algorithm can be trivially extended to optimize (both to minimize or 
maximize) other parameters of TZ'Hg{P), such as the perimeter (the sum of the length of 
the boundaries of all connected components) or the number of connected components of 
TPHq{P). Interesting extensions of the “parameter optimization” family of problems, are 
to consider dynamic or moving point sets. In particular, the characterization of the area 
function on each variation seems an interesting and non-trivial problem. 
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